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The MINOS long-baseline experiment is using the NuMI neutrino beam to make pre- 
cise measurements of neutrino flavor oscillations in the "atmospheric" neutrino sector. 
MINOS observes the disappearance oscillations seen in atmospheric neutrinos, tests 
possible disappearance to sterile v by measuring the neutral current flux, and extends 
our reach towards the so far unseen 613 by looking for appearance in this Vf^ beam. 
The magnetized MINOS detectors also allow tests of CPT conservation by discrimi- 
nating between neutrinos and anti-neutrinos on an event-by-event basis. The intense, 
well-understood NuMI neutrino beam created at Fermilab is observed 735km away at 
the Soudan Mine in Northeast Minnesota. High-statistics studies of the neutrino inter- 
actions themselves and the cosmic rays seen by the MINOS detectors have also been 
made. MINOS started taking beam data in May of 2005 and is now nearing the end of 
it's five-year run. This paper reviews results published based on the first several years of 
data. 
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1. Introduction 

One possible implication of non-zero neutrino mass is that as neutrinos propa- 
gate, mixing occurs between mass eigenstates (1^1, i^2, 1^3) and the flavor eigenstates 
[vej u^j^^Vr) involved in the weak interactions governing neutrino production and de- 
tection. The mass and flavor eigenstates are related by a unitary matrix Upmns 
expressed in terms of three mixing angles (6*12,^13,023) and a CP-violating phase 
5 [l]l2] . Experimental evidence points towards this occurring via disappearance of Ve 
in solar neutrinos [3l|4] and D^. in reactor neutrinos [5], as well as i'^ disappearance 
in atmospheric neutrinos produced by cosmic rays [BHS] and accelerators |9]. 

The goal of the Main Injector Neutrino Oscillation Search ("MINOS") long- 
baseline experiment is to precisely measure the oscillation parameters involved in 
the atmospheric-sector oscillations. It does this by observing the intense and well- 
understood NuMI beam over a known baseline using with two similar magnetized 
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steel/scintillator calorimeters: 1 km from its source at Fermilab with a 0.98 kton 
"near detector", then again 735km to the northwest in the Soudan Mine Under- 
ground Lab using the 5.4 kton "far detector". Such a before-and-after comparison 
of the neutrinos greatly reduces systematic errors associated with comparing dif- 
ferences in the observed neutrino spectra to various neutrino oscillation scenarios, 
allowing for a more accurate probe of the physics of neutrino propagation. Details 
of the detector designs, calibrations, and performance can be found in [lOj . 

The NuMI beam [llj is generated by 120 GeV protons hitting a carbon tar- 
get. The resulting charged pions are focused by two electromagnetic horns and 
sent down a 675 m decay pipe, producing a beam of 92.9% i/^, 5.8% F^, 1.2% i^e 
and 0.1% Fg. Changing the pion-focusing horn positions and currents creates very 
different neutrino spectra. The bulk of the data come from the "low energy" beam 
configuration, peaked at several GeV (see the dashed line in Fig.[T]). Short exposures 
in other beam configurations observed by the high-statistics near detector provide 
good crosschecks to the beam modeling process and also help to reduce systematic 
errors. 

This paper summarizes the results of several analyses of the neutrino data ac- 
quired over the two year time period starting with the beginning of NuMI operations 
in May of 2005 and ending during the summer shutdown in June 2007, an integrated 
exposure of over 3 x lO^'^ protons on target ("pot") with a neutrino yield on order of 
one neutrino per proton. The intrinsic divergence of the beam results in a neutrino 
flux at the far detector which is a factor of 10^ lower than that at the near detector. 
An additional 4 x lO^'' pot is on tape, representing beam from October 2007 through 
June 2009. The beam focusing was reversed to favor anti- neutrino production upon 
resumption of data taking in September of 2009, and was switched back to normal 
neutrino mode in March 2010 after collecting 1.8 x 10^" pot. 

MINOS analyses are done "blind" . That is, small subsets of the data are ex- 
amined immediately to develop the algorithms and cuts which extract the physics 
results from the data, but the bulk of the dataset is left alone till the data analysis 
methods have been settled on. Detailed comparisons of variables used in cuts are 
made between simulated Monte Carlo data and real data, to be sure that detec- 
tor quirks and parameter distributions are well modeled and understood. Once the 
simulations and analysis techniques are close to final and uncertainties, systematic 
errors and expected sensitivities calculated, "sidebands" (parameters similar to but 
not actually the ones used to produce the final answer) are revealed and compared 
to expectations. Only once everything is as well understood as possible is the "box 
opened" and the mature analysis unleashed on the full dataset to produce the final 
physics results. This process helps to ensure that the development of the analysis 
doesn't inadvertently (consciously or otherwise) converge on a fluctuation in the 
data, skewing the flnal results or their significances. Extensions to the analyses 
presented below to the newer, as yet unanalyzed data are in progress, and as new 
improvements are being made to the methods, the new analyses are also using this 
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blinding procedure. 

2. Oscillation Analyses 

MINOS was intended to observe Charged Current ("CC") quasi-elastic i/^ inter- 
actions to make a precision measurement of Amgj. Its calorimeters are designed 
to measure the fi produced in CC interactions, and are magnetized to provide 
momentum and charge discrimination on an event by event basis, as well as the 
standard calorimetric method of establishing particle energies via dE/dx. By using 
far and near detectors which are as similar as possible, uncertainties in the energy 
scale between measurements of the neutrino spectrum before and after traveling 
the 735 km baseline are minimized. Small differences between detectors in light col- 
lection and electronics were cross-calibrated in a beam test at CERN using the 
"calibration detector" [T^] which observed the same particles with both sets of elec- 
tronics P^ . 

Since the oscillation minima at this baseline is less than the t production thresh- 
old energy, the oscillatory signature is that of j/^ disappearance rather than Vr ap- 
pearance, an analysis described in Sec. 12.11 MINOS' magnetic fields allow discrimi- 
nation between and and provide a chance to test CPT conservation by seeing 
if neutrinos and anti- neutrinos have the same oscillation parameters fSec. 12.2"]) . Al- 
though the calorimetry is coarser than one might like for such a measurement, 
electromagnetic showers can still be resolved. This allows further study of possible 

disappearance to sterile ly states by observation of Neutral Current ( "NC" ) inter- 
actions fSec. [273)) . and sensitivity to sub-dominant i/^ o transitions by looking 
for appearance oscillations (Sec l2.4p . 

2.1. Ufj, Disappearance Oscillations 

Since 623 3> ^13 and 612, when considering the atmospheric neutrino sector the full 
neutrino mixing matrices reduce to a simpler two-flavor equation. A of energy 
£'i^[GeV] observed after traveling some distance L[km] from its production point 
has a probability of being detected as a i^f^ given by 

P(j/^ -^u^) = l- sin2(2(?23) sin2 (^1.27Am22|^ , (1) 

where Am^ [eT^^] is the mass difference between i>2 and and sin^ (26) is the mixing 
amplitude. 

An exposure of 3.36 x 10^" pot in MINOS has been analyzed [H], selecting 848 
far detector events as with good purity. The observed (unoscillated) near detec- 
tor signal is used to calculate a null hypothesis expectation of 1065±60 far detector 
events, including a small background estimated to be composed of 2.3 external /x, 
5.9 NC induced showers, and 1.5 r decays. The resulting spectra is shown with the 
observed data in Fig. [T] Comparing the observed data to expectations modified by 
Eq.[T] result in best fit oscillation parameters of |Ato^| = (2.43 ± 0.13) x 10~^ eV^ 
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Fig. 1. The MINOS far detector spectrum 14 . The data (points with statistical errors) show 
a significant deficit from the null hypothesis (dashed line), but well- match a -fl oscillation 
scenario (solid line), with best fit mass splitting |Am2| = (2.43 ±0.13) X 10"^ eV^ (68% cl) and 
mixing angle sin2(2e) > 0.90 (90% cl). 

(68% cl) and sm'^{2e) > 0.90 (90% cl). Systematic errors in Am are dominated 
by uncertainties in the hadronic energy scale (±10.3% absolute and ±3.3% relative 
between near and far detectors) and the relative normalization between detectors 
(±4%). The background of NC showers mis- reconstructed as quasi-elastic in- 
teractions (±50%) dominates in sin^ 2923- These systematics (±0.108 x lO^^eV^ in 
Am^ and ±0.018 in sin^ 2823) are still smaller than the statistical (±0.19 x 10"'^ eV^ 
in Am^ and ±0.09 in sin^ 2^23) errors, so the measurement will improve as data 
from the remaining ^ | of the MINOS exposure are added. 

2.2. Anti-neutrinos 

The magnetized nature of the MINOS detectors allows the event- by- event determi- 
nation of the charge sign of muons, and thus the identification of the parent neutrino 
or anti-neutrino undergoing the quasi-elastic interaction that produced them. Se- 
lection of wrong-sign muons in the beam tests if t'^, oscillate in the same fashion 
as in Eq. [H Is 623 = O23 and Am§2 = Arn|2? Furthermore, could some fraction 
a of disappearing j/^ reappear as These are both tests of CPT conservation. 
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Fig. 2. The resulting allowed region in oscillation parameter space for the data in Fig. [T] at 90% 
(solid black) and 68% (dashed black) confidence levels. Two Super-K atmospheric analyses are 
shown in red [7] and blue [6], and the K2K long baseline experiment results produce the gray 
contour [9]. 



and with its magnetic field MINOS is the only experiment capable of testing this. 

Anti-neutrinos in the NuMI beam come primarily from n'^ which travel directly 
down the center of the focusing horns and thus avoid being de-focused, leaving little 
kinematic phase space for these pions. Combined with the lower cross-sections for 
anti-neutrinos compared to neutrinos, only 6.4% of the neutrino interactions in a 
3.2 X 10^° pot far detector exposure are due to anti-neutrinos, making the relative 
backgrounds are higher, the statistics lower, and the spectrum harder (Fig. [3]). 
42 events are seen while 64.6 ± S.Ostat ± S.Osy^t are expected in the no-oscillation 
case, or 58.3 ± T.Qstat ± 'i-&syst if CPT is conserved given the observed oscillation 
parameters [15]. This places a 90%cl upper limit of a < 0.026, and the anti-neutrino 
oscillation parameters are consistent with the neutrino parameters given these low 
statistics, excluding 5.0 < Amgj < 81 eV^ at 90%cl for maximal mixing. 
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Fig. 3. The MINOS far detector spectrum [15J. The data (points with statistical errors) are 
consistent with CPT-conserving neutrino oscillations (blue dashed line) within the large error bars. 

To increase the anti-neutrino dataset by an order of magnitude, the polarity 
of the focusing horns in the NuMI beam was reversed to focus 7r+ rather than the 
usual 7r~ for an exposure of 1.8 x lO^'^ pot of anti-neutrino production. The statistics 
available in this set of data would make a reasonable measurement of 623 and Am^2 
(FigHl). 



2.3. Sterile Neutrinos 

Another possible explanation of i'^ disappearance is oscillation into sterile neutrinos 
("fs") which experience no interactions and thus would disappear. This would also 
suppress the rate of NC events in the far detector compared to the traditional expla- 
nation of oscillation to sub-threshold Vr, since those h'r still undergo NC interactions 
but would not. 

To test this hypothesis, NC showers have been selected from an exposure of 
3.18 X lO^'^ pot [H]. A NC interaction produces no outgoing leptons, but simply a 
hadronic shower. Any tt" produced in that shower decays rapidly to a pair of gamma 
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Fig. 4. Sensitivity contours comparing the oscillation parameters available using just the anti- 
neutrinos from the neutrino beam, as described in the text but including all available data 
(7.2 X 10'^'^ pot, the black line) to that obtainable using a 7.2 X 10^" pot exposure of dedicated 
anti-neutrino beam (red contour). This red contour is comparable in scope to that of the K2K 
experiment's neutrino results (the gray line in Fig.[2]|. 

rays, which produce diffuse electromagnetic showers that can be reliably separated 
from the long /i tracks used in Sec. 12.11 388 events are selected, 141 of them of 
less than 3GeV, with an estimated 17 non-NC interactions (primarily very short 
track CC interactions) creeping in as background in this low energy region of 
interest for oscillation physics. The resulting spectrum of these NC events (Fig. [5|) 
is not depressed compared with the expectations of NC interactions from standard 
i^fj. ^ i^T oscillations. The ratio of observed to expected NC events in the far detector 
is i? = 0.99±0.09s(at±0.07syst — 0.08^^ in the energy region of 0-3 GeV where are 
disappearing. While the addition of a fourth neutrino eigenstate results in several 
new ways things could oscillate, fits to all those models place limits on 024 and 6*34 
on being half (or less) than 623, showing that sterile neutrinos are not a dominant 
player in disappearance. Another way to state this is as a limit on the fraction 
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of Vs participation of fs < 0.51 at 90%cl. 



4 — Far Detector Data 




Ereco (GeV) 

Fig. 5. The MINOS far detector neutral current spectrum |16| . The data (points with statistical 
errors) track the predicted spectrum (red hashed boxes) and do not show the deficit at lower 
energies expected from an scenario which involves the NuMI changing to something which does 
not undergo such interactions, either a sterile neutrino or by decaying away entirely. The dashed 
line is the prediction for showers should there be Ug appearance at the CHOOZ limit, as they would 
be included in this data selection. The gray shaded region is the expected background, primarily 
composed of mis-reconstructed charged current events. 



2.4. The Search for Uf. Appearance 

The CHOOZ experiment [17] sets an upper limit of sin^(26'i3) < 0.15 on the mixing 
amplitude governing the transmutation of NuMI into v,,. MINOS was designed 
to be a good muon calorimeter for disappearance, but is coarse for resolution of 
~GeV electromagnetic showers, each 2.54cm steel+l.O cm plastic thick plane being 
1.4 radiation lengths thick, and each 4.1 cm wide scintillator strip being 1.1 Moliere 
radii. However, the experiment retains sensitivity to the ~2% v,, appearance signal 
which a ^13 near the CHOOZ limit [T7| would create, and the first 3.14 x lO^Opot of 
MINOS data have been examined [18] by a neural network to select electromagnetic 
shower candidates, which are more compact when produced by an electron from a 
CC Ve interaction than from a NC-induced 7r°. When applied to Monte Carlo data 
this is 41% efficient at keeping Ue events while rejecting >92% of NC showers (the 
dominant background) and >99% oiv^ charged current ("CC") interactions (high-y 
collisions which put much more energy into the hadronic debris than the outgoing 
(1 lepton). 



April 16, 2010 0:8 WSPC/INSTRUCTION FILE minos 



MINOS neutrino oscillation results 9 

Given the small expected signal and large uncertainties in hadronic shower mod- 
eling, data-driven methods are used to better estimate the background. At the near 
detector no oscillation has yet occurred, so with the exception of the well-modeled 
inherent beam Ve, all events selected must be examples of such background events. 
The two classes of backgrounds have different production kinematics so extrapolate 
to the Far Detector slightly differently, thus two techniques are used to deconvolve 
the background. The first takes obvious interactions and subtracts the hits from 
the muon track, resulting in a sample of CC-induced showering events to study. The 
second compares data from beam running with the focusing horn on or off. The very 
different neutrino spectra which result allow fitting for the two background compo- 
nents. Both methods give comparable results and produce the background at the 
far detector shown as the red line in Fig. [6l This yields an expected background 
of 26.6 (18.2 NC, 5.1 CC, and 2.2 beam I'e) at the far detector, while 35 like 
events are seen, a 1.5 a excess (including 7.3% statistical and 19% systematic errors) 
(Fig. [6|). If fit for oscillations, this is just below the CHOOZ limit and consistent 
within errors with no Ve appearance. 
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Fig. 6. The spectrum of potential Ue interactions in the MINOS far detector, with statistical plus 
systematic error bars |18| . The l.Scr excess is consistent with both the expected large background 
(red) and a sin^(29i3) comparable to the CHOOZ limit |17| (purple). 



There is somewhat more than twice the exposure used for this analysis already 
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on tape. In addition to decreasing the statistical error on the data points, the 
systematic error bars include a large statistical component due to the small data 
subsets used in the data-driven background estimation. Thus, an analysis of the 
additional data available is expected to substantially reduce the systematic error 
bars as well. Projections of the sensitivity of this new dataset suggest that the slight 
excess observed in this current analysis will either be shown to be more significant 
(if it is really Ve appearance) or to be revealed as merely a background fluctuation. 

2.5. Alternate Hypotheses 

Two other hypotheses have been presented (and not thoroughly ruled out by previ- 
ous observations) to explain the disappearance outside the standard oscillations 
model. The first is quantum decoherence of the neutrino's wave packet [20[ . in which 
the survival probability equivalent to Eq.[T] varies as [1 — exp(2^^^)] (Eq. 5 of [20]) 
rather than sin'^ {1 .27 L / E) . The best fit to this function is shown in Fig. [7] as the 
dashed line, and a comparison of to the data and the standard O i^r disfavors 
the decoherence hypothesis at the 5.7 a level. 

The second alternate hypothesis is that neutrinos decay [19 , producing a sur- 
vival probability which varies as [sin^ 9 + cos^ 6'exp(^^^)]^ (Eq. 13 of [19 ). Again, 
this shape does not match the data in Fig. [7] as well as standard oscillations, and 
is disfavored at 3.7 a. However, this result is for pure neutrino decay, which is also 
disfavored by Super-K [6]. If there were oscillations and neutrino decay happening 
at the same time [21], the waters are muddier, see Eq. 17 of [TIT. In order to get 
a better handle on this problem, MINOS makes use of NC events as well as the 
CC spectral shape [16] , since a decayed neutrino will not produce a NC event while 
a Vt will. Doing this improves the rejection of the pure decay hypothesis to 5.4 u, 
and places a limit on the ratio of neutrino mass eigenstate rria's lifetime to mass of 
ra/mg > 2.1 x lO^^^s/eV (Fig. [8]). 

2.6. Atmospheric Neutrinos 

In addition to the NuMI beam, the MINOS far detector is bathed in the same 
fiux of cosmic-ray induced atmospheric neutrinos which provided the first measure- 
ments of neutrino oscillations [22 . However, MINOS being an order of magnitude 
smaller than Super-Kamiokande, the atmospheric neutrino interaction rate is sev- 
eral per week rather than dozens per day, limiting the statistical significance of such 
measurements. Nevertheless, analyses of such neutrinos are consistent with the os- 
cillation parameters established in the beam neutrinos [23] and provide the first 
direct observations of anti-neutrinos from cosmic rays [33] . 

3. Non-Oscillation analyses 

This review is of MINOS' neutrino oscillation results, so other work using these 
detectors will not be discussed in depth. However, studies using cosmic rays have 
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Fig. 7. The ratio of the MINOS far detector data from Fig.[T]to the no-oscillation null hypothesis, 
from |14| . Superimposed are the best fit expectations of three disappearance models: standard 
Vp^ •<-> Ur (thick solid line), which fits the data the best; pure neutrino decay |19l (thin solid 
line), which is disfavored at the 3.7cr level; and quantum decoherence [20) (dashed line), which is 
disfavored by this data at 5.7 cr. 



grown out of the need to calibrate the MINOS detectors, and speak to the depth 
of understanding of these detectors. Cosmic ray analyses include the first direct 
measurement of the charge ratio of cosmic ray muons at TeV energies [25j and 
probe the meson production in cosmic ray primary interactions by watching the 
variation in the underground muon rate vary with stratospheric conditions |26j . a 
competition between secondary mesons decaying to produce the observed cosmic ray 
muon and re-interacting in the atmosphere. This effect is also of use to atmospheric 
physicists, who turn the problem around to use the cosmic rays to study unusual 
events in the stratosphere itself [27^. 

Studies of neutrino interactions themselves are also a large topic of work, as 
the near detector observes 0(10'^) neutrino interactions per day of operations, by 
far the largest statistics sample in the world. The resulting improved knowledge of 
neutrino interaction physics feeds directly back to reducing the systematic errors 
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Fig. 8. MINOS' 90%cl allowed region for the neutrino mass/lifetime ration a compared to the 
oscillation mixing angle 9 (T6| . The best-fit is for no decay (infinite lifetime), with an upper limit 
of T3/m3 > 2.1 X 10-12 g/gv. 

in the oscillation analyses p8|. Looking at MiniBOONE neutrinos in MINOS and 
vice- versa helps to understand the ofF-axis components of neutrino beams [29j. Using 
beam neutrinos in different ways has also yielded interesting results on the velocity 
of neutrinos (comparing arrival times at near and far detectors) [30] and has been 
used to test for violation of Lorentz Invariance in the neutrino sector [31j . 

4. Conclusions 

MINOS has measured neutrino oscillation parameters in the "atmospheric" 1^2 ^ 
sector with high precision, favoring standard t/^ -s-> Vt oscillations with |Am^| = 
(2.43 ± 0.13) X 10-3 eV2 (68% cl) and sin^{2e) > 0.90 (90% cl). Quantum deco- 
herence as an explanation for the disappearance is disfavored at the 5.7 a level. 
Measurements of the total active neutrino flux using neutral current interactions 
help to disfavor pure ^ Vs oscillations by 5.4 ct, place a limit on the ratio of neu- 
trino mass eigenstate ma's lifetime to mass of T^/mj, > 2.1 x IQ-^^s/eV, and limit 
the participation of sterile neutrinos in a sub-dominant mode to a fraction fs < 0.51 
at 90%cl. Examination of the limited set of anti-neutrino data have shown no evi- 
dence for different oscillation parameters for anti-neutrinos, and limit the fraction 
a of Vfj, disappearing to to a < 0.026 at 90%cl. 

These results come from an exposure of roughly 1 /3 the eventual complete MI- 
NOS dataset, so the final precision of the measurements will be greater than those 
reviewed here. This is especially true in the cases of appearance, where the 
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data-driven background estimations methods will benefit greatly from additional 
statistics, and in the measurements of P^, which will take advantage of dedicated 
anti-neutrino beam running. 
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